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Abstract. We implement the dressing method for a novel integrable gen- 
eralization of the nonlinear Schrodinger equation. As an application, explicit 
formulas for the A^-soliton solutions are derived. As a by-product of the analy- 
sis, we find a simplification of the formulas for the A^-solitons of the derivative 
nonlinear Schrodinger equation given by Huang and Chen. 
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1. Introduction 

We consider the following integrable generalization of the nonlinear Schrodinger 
(NLS) equation, which was first derived in [2] by means of bi-Hamiltonian meth- 
ods: 



> 



o 



^1.1) iu t — vutx + iu xx + a\u\ (u + iuu x ) = 0, a = ±1 



where 7 and v are nonzero real parameters and u(x,t) is a complex-valued 
function. Equation ( II .ip reduces to the NLS equation when v — 0. However, 
the limit v — > is in many respects singular, and the analysis of (II .ip is rather 
■ different from that of NLS. Equation (II .ip appears as a model for nonlinear pulse 

! propagation in optical fibers provided that one retains certain terms of the next 

asymptotic order beyond those necessary for the NLS equation (2j [5]. From a 
bi-Hamiltonian point of view, equation ( II. ip is related to NLS in the same way 
that the celebrated Camassa-Holm equation [T] is related to KdV [2J. Being 
integrable, (II. ip admits a Lax pair formulation and the initial- value problem on 
the line can be analyzed by means of the inverse scattering transform (1ST) [6]. 
The spectral analysis of (11.11) is closely related to that of the derivative NLS 
equation. In fact, equation (11.11) is related by a gauge transformation to the 
first negative member of the integrable hierarchy associated with the derivative 
NLS equation [5] . The initial-boundary value problem for equation (II. ip on the 
half-line was studied in [TJ. 

Here, we implement the dressing method for equation (11.11) and derive, as 
an application of the general approach, explicit formulas for the iV-soliton so- 
lutions. The dressing method is a technique which generates new solutions of 
an integrable equation from an already known seed solution cf. [8]. In our 
implementation of the dressing method to equation (ll.ip . we start with a seed 
solution uq together with a corresponding eigenfunction ipo of the Lax pair. A 
new eigenfunction ip = Gipo is constructed by multiplying ip by a 2 x 2-matrix 
G which contains an even number of poles in the spectral parameter. It is shown 
that ip satisfies the same Lax pair as ipo but with the potential Uq replaced with 
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a new potential u(x, t). The compatibility of the Lax pair equations then implies 
that u(x,t) also is a solution of (II. ip . 

An expression for the one-soliton solution of ( 11. ip was found in [6] by solving 
directly the Riemann-Hilbert problem associated with the inverse problem in 
the presence of two poles and no jump (each soliton contributes two poles to 
the Riemann-Hilbert problem because the poles are of necessity symmetrically 
distributed with respect to the origin). Applying the dressing approach to the 
particular seed solution u = we recover the one-soliton solution and also find 
an explicit expression for the A-soliton solution for any N. 

A convenient Lax pair for (II. ip is introduced in Section [2j while the implemen- 
tation of the dressing method is described in detail in Section [3j In Section H] we 
derive explicit expressions for the iV-solitons. Finally, in Section [5j motivated 
by the analysis of (11.1 p . we find a simplification of the formula for the iV-soliton 
solution given by Huang and Chen [3] of the derivative nonlinear Schrodinger 
(DNLS) equation 

(1.2) iqt + qxx + i(\q\ 2 q) x = 0, x e R, t > 0. 

The connection between equations (II. ip and (II. 2p arises because the x-parts of 
the corresponding Lax pairs are identical upon identification of u x with q. 



2. A Lax pair 

It was noted in [5] that the nonzero values of the parameters 7, 1/, a in (II. ip 
can be arbitrarily assigned via a change of variables. We will therefore, for 
simplicity, henceforth assume that 7 = v — a — 1. Then the transformation 
u e lx u converts (II. ip into the equation 

(2.1) Utx + u — 2iu x — u xx — i\u\ 2 u x = 0, 

which admits the following Lax pair [6]: 

{ip x + i-ko- 3 ip = tU x iP, 
2 
A + i - |) ^ = (\U X - \o z U 2 + f^tf) V, 

where ip(x,t,() is a 2 x 2-matrix valued eigenfunction, ( £ C = CU is a 
spectral parameter, and 

We will assume that ip obeys the symmetries 

(2.4) a 3 il>(x,t,-Oa 3 =i/>(x,t,C), ip- 1 {x,t,() = ^{x,t,(). 

A solution ip of (I2.2p satisfying (12.41) can be constructed as follows: Assuming 
that u has sufficient smoothness and decay, there exists a solution of (12.21) which 
tends to exp(— i6cr 3 ) at infinity, where 

(2.5) 0:=e(x,t,C) = ^x+ Q-0 t. 
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This suggests the transformation if) = \l/e %9(T3 and then \1/ solves 
J * a , + i^r[<T 3 ,*] = ±U X V, 

(2 ' 6) | * t + i(i-|) 2 [<7 3 ,*]= (iC/ s -itT 3 C/ 2 + |(73C/)*. 

We define two solutions and ^2 of (12. 6ft via the linear Volterra integral 
equations 

(2.7a) tt 1 (z,t,C) = /+ f e^^-^C"^^*!)^',*^)^', 

«/ —00 

/■oo 

(2.7b) * 2 (x,t,C) = J- / e^'-^C^^X^Oeb', 



a; 



where <r 3 acts on a 2 x 2 matrix A by o" 3 y4 = [(73, A]. The second columns of 
these equations involves exp[2i£~ 2 (x' — x)]. It follows that the second column 
vectors of \&i and \I/ 2 are well-defined and analytic for ( G C such that £ 2 lies in 
the upper and lower half-planes, respectively. Moreover, these column vectors 
have continuous extensions to £ G R U iR \ {0, 00}. Similar remarks apply to 
the first column vectors. 

We define a solution if) of ( 12. 2ft by 



(2.8) ^(M,C) 



([* 2 ]i, [*i] 2 )e- i0<T3 , ImC 2 >0, 
([*i]i,^2] 2 )e-^ 3 , ImC 2 <0, 



where ([ x f , 2 ]i, [ x ^i]2) denotes the matrix consisting of the first column of \I> 2 
together with the second column of \&i etc. Then if) admits the symmetries in 
(12.41) . Indeed, uniqueness of solution of (I2.2p implies that there exist matrices 
Si(C) and S 2 (C) independent of x, t such that 

a 3 ifj{x, t, -C)a 3 = ifj{x, t, C)Si(C), ^{x, t, () = S 2 {()^{x, t, £), 

and evaluation of these equations as x — > ±00 using (12. 7p shows that Si = S 2 = 
I. 

As a function of (, the eigenfunction if) defined by (12.81) is singular for £ e 
{0, 00} and has a jump across the contour R U iR of the form -0_ = ip + J(£), 
where ip± denote the values of if> on the left and right sides of the contour and 
J(C) is a 2 x 2 jump matrix independent of x and t. However, these singularities 
are inconsequential for the arguments below (which use the combinations if) x if)~ l 
and iptip' 1 )- I n fact, Liouville's theorem implies that if>, in general, cannot be 
analytic over the whole Riemann (^-sphere. 

3. The dressing method 

Starting from a seed solution uq of (12. ip and a corresponding eigenfunction ip 
obeying the symmetries ( 12.41) f*l we seek a 'dressed' eigenfunction if) of the form 

o- 



(3.1) ip = Gip, 



1 As outlined in the previous section, such a ipo can be constructed from uq using only linear 
operations. 
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where the 2 x 2-matrix valued function G has the form 

N r 

(3.2) G(x,t,C) = I + Yl 



Aj(x,t) (T z Aj{x,t)(T 3 



c-o 



C + Cs 



Here {Cj, — C?}i CC\{0, 00} is a collection of simple poles with corresponding 
residues {Aj(x, t), —a 3 Aj(x, t)a 3 }^ . The form (13. 2p of G is motivated by the 
condition that G — > I as ( — > 00 and by the symmetry 

(3.3) a 3 G(x,t,-C)cr 3 = G(x,t,0, 

which ascertains that the first symmetry in (I2.4p is preserved by (13.11) . Moreover, 
in order for the second symmetry in (12. 4p to be preserved by (13. ip . we require 
that 

(3.4) G- 1 (x 1 tX) = G^x } t } () } 



i.e. 



(3.5) 



N 



G- 1 (x,tX) = I + Yl 



A}(x,t) a 3 A){x,t)a 3 



c-o c-o 

It follows from (13. 2p and (13. 5p that G and G~ l are analytic at both Q = and 

C = 00. 

The goal now is to show that ip satisfies (12 .2p for some matrix U of the form 
(12. 3p . The compatibility of the Lax pair equations will then show that u := U12 
is a solution of equation (12. ip . 

Differentiation of (13. ip with respect to x and t followed by multiplication by 
= ^p 1 G^ 1 leads to 

(3.6a) = G X G- 1 + G^o'G' 1 

and 

(3.6b) ^VT 1 = G t G~ x + G^o'G- 1 , 

respectively. 

3.1. Analysis at £ = 0. We replace ^Oxipo 1 and ipotipQ 1 in (13 .6p with the fol- 
lowing expressions which follow from the Lax pair equations satisfied by ipo'- 



(3.7) 



2 



Mo 1 = ~i (j - 3)^3 + \U 0x - ia 3 Ul + f^C/o- 
Substituting into the resulting equations the expansions 



G = Go + GtC + 0(C), G~ = G L - G^dGoX + Q(( 



where Go(x,t) and G\(x,t) are independent of (, it follows that ipxip -1 and 
iptip -1 have double poles at ( = 0; identification of terms of 0(l/( 2 ) and 0(l/£) 
yields 



(3.8) 



^- 1 = % + % i + o(i), 



C 2 ' c 

Q-2 , g_i 



c 5 



c 



0(1), 
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where 

(3.9) g_ 2 = -iG Q a 3 Gu\ 

Q-i — iGocr^GQ G X G 1 — iGi<j 3 G 1 + GqUq x G . 
On the other hand, from ( 13 .2p we find the following expressions for Go and G\\ 

n _ r A 3 + gV^Tg _ ^ - OsAjOs 

G o-^-2^ 7 ' ^-"Z^ 72 ' 

J=l ^ i=l ^' 

so that the matrices Go and G\ are diagonal and off-diagonal, respectively. Thus 
the expressions for Q_ 2 and Q-i i n (I3.9j) simplify to 

(3.10) Q_ 2 = -itra, Q-i = -2iG x Gz x <Ts + GoC/o.Cq- 1 . 
In particular, Q_i is an off-diagonal matrix. 

3.2. Analysis at £ = oo. We again replace ipox^® 1 and tpotipQ 1 in (I3.6P with the 
expressions in (13. 7p . Substituting into the resulting equations the expansions 

G = / + ^ + ^ + 0(l/C 3 ), C^oo, 

where G_i and G_ 2 are independent of (, it follows that ipxip' 1 — >■ as £ — > oo 
whereas ipttp -1 has a double pole at £ = oo; identification of terms of 0(( n ), 
n = 0, 1, 2, yields 



(3-11) G- 1 = /-^ + "- 1 ,, ~ 2 +0(l/C 3 ), C^oo, 



= ->sC 2 + QiC + Qo + 0(i/C), C oo, 



where 



+ ia 3 - u -(JzUl. 



(3.12) Q 1 = l[a 3 ,GU] + ^X3f/o, 

11% 1 1 

Qo = -^W3,G- 2 ] - -jO-&G\ + -G_ia 3 G_i + -G_ier 3 f7 - -^^oG'-i 

z 

— t 
2 

On the other hand, from ( 13. 2 p we find the following expressions for G_i and 
G_ 2 : 

2V AT 

(3.13) GU = ]T (^ - a 3 A j( T 3 ) , GL 2 = ^ Q (A, + osA,-^) , 

i=i i=i 

so that the matrices G_i and G_ 2 are off-diagonal and diagonal, respectively. 
Thus the expressions for Qi and Qo in ( 13. 12ft simplify to 

(3.14) Q x = l -a,U, Qo = ~U 2 a 3 + ia 3 , 
where we have defined 

(3.15) U := GLi + C/ . 
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Equations (I3.5P and (13. lip yield the following alternative expression for GLi: 

N 



Thus (GLi)2i = — (G_i) 12 , so that U is of the form (12.31) with u := (G-i)i2 + Mo- 



3.3. The generated solution. It follows from the previous two subsections 
that the functions 

(3.16a) M- 1 - ( % + Q - 



e C 



and 



(3.16b) w 1 -l^ + ^f + Qo + QiC- 

are analytic near ( = and £ = 00. Thus, in view of (13. 6p and the expressions 
(13. 2p and (13. 5p for G and they are analytic on the whole Riemann ^-sphere 
except possibly at points in the set {±Cj> ^Cj}i ■ We claim that the A/s in (13. 2p 
can be chosen so that the functions in (13.161) are analytic everywhere. Assume 
that such a choice has been made. Then, since the functions in ( 13 . 1 6[) tend to 
zero as ( — >■ 00, Liouville's theorem implies that they both vanish identically. 
Together with the expressions (I3.10p and (13. 14j) for the coefficients Q-2,Qo,Qi, 
this implies that ip satisfies the following pair of equations: 

2 

& + i (? - §) a ^ = (jQ-i - \^U 2 + fau) ^ 



where U has the form (I2.3p . Using that Q-\ is off-diagonal and identifying terms 
of 0(C) m the compatibility equation ip xt = ip tX) we find that 

(3.18) Q_! = U x . 

Consequently, (I3.17P has exactly the form of the Lax pair (12. 2p and its compati- 
bility implies that the function u = U12 satisfies (12. ip . We infer from (I3.13P and 
(I3.15P that the generated solution u(x,t) is given by 

N 

(3.19) u(x,t) = 2^(A j (x,t)) 12 + u {x,t). 



3.4. The dressing transformation. In order to complete the construction of 
the map Uq h-> u, it only remains to determine the structure of the A/s in (I3.19P 
so that the functions ip x ip~ l and iptip -1 are regular at the points 
To this end we consider the following series of N consecutive dressing transfor- 
mations, each of which adds two poles: 



(3.20) G = D N D N -! ■■■D 1 , ^ = ^-1, j = l,...,N, 
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where, for j — 1, . . . , N, 
(3.21a) Dj(x,t,C) = 1+ 3 



(3.21b) Dj\x,t,()=I + 



c-o c+o ' 

B]{x,t) a 3 B](x,t)a 3 



c-o C+0 

and {-Bj(x, are a set of 2 x 2-matrix valued functions independent of (. 

Lemma 3.1. Let {bj}i be a set of N nonzero complex constants. Define the 
functions Bj(x,t), j = 1, . . . , N, inductively by 

(3.22) Bj(x,t) = \ Zj {x,t)}{ yj {x,t)\, 

where the column vector \zj) = (zj\* and the row vector (yj\ = \yj)^ are defined 
in terms of the (j — l)th eigenfunction ipj-i by^ 

(3.23) (y j \=(b j bj 1 )^-^), 

(3-24) : ; ) (']; !!) // ; ). where a? = < % | ") <;,). 



Define ip := ip^ by li3.20\) - (3.24\ ). Then the functions ip x ip~ l and iptip -1 are 
analytic at the points in the set {±Cj> ^Cj}i ■ 

Proof. The proof proceeds by induction. Suppose that has been defined by 
f l3.20p - fl3.24p and that the functions ipj-i^ipj^ and ipj-ijipj^ are analytic at 
( = ±Q and ( = ±Q for i = 1, . . . ,j — 1. We will show that and ipjtipj 1 , 

where ipj = Djipj_i, are analytic at £ = ±Q and ( = ±Q for i = 1, . . . , j. 
Differentiation of ipj = Djipj_i yields 

(3.25) V;^" 1 = D jx Dj l + I),c, \D] '. 

The assumption on ipj-i together with ( 13.211) implies that the right-hand side is 
analytic at the points {±C*j ±Qy{~ . To show that ipjxipj 1 is regular at Q, we 
note that ( 13. 25ft implies that ipjxipj 1 has at most a simple pole at Q, and that 

(3.26) Res (<.>,.•') = {B^^Q^j^QDj 1 ^). 

Using fl3T22|) and fl3T23|) . we can write fl3T26|) as 

Res (ipjxipj 1 ) = \z j ) x {y j \Dj 1 {Q. 

We claim that 

(3-27) (Vi\Dj%) = 0, 

so that the residue vanishes. Indeed, by ( 13.2 lbj) . 

(3.28) MDF %) = fc| + - <fJfMfeK 

and the right-hand side of this equation vanishes by virtue of f)3.24p . Since 
the symmetry properties fl2.4[) are preserved by (13.201) . we deduce that ipj X ipj 



2 Hcrc and in some equations below the (x, ^-dependence is suppressed. 
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is analytic also at — Q and ±(j. Similar arguments establish the regularity of 

VjtVj 1 . ' " □ 

By (USD and (13T201 . 

(3.29) A j = ResG(C) = D N (Q . . . D^QB^^Q . ..D 1 (Q). 

Q 

On the other hand, it follows from Lemma 13.11 that 

(3.30) a 2 Bj(x,t)a 2 = C -t^Dr\ x ,t,Q, * 2 = ~fj , 
and 

c 2 - c 2 

(3.31) a 2 D 3 (x,t,() T a 2 = l^-Dj^xM)- 

Using (I3.30p and (13.311) . the expression (I3.29|) for Aj can be written as 

(3.32) AAx, t) = -a 2 G-\x, t, Q T a 2 



a i 



where the complex constants {ctj}i are defined by 
Substitution of (I3.32j) into the expression (13. 2p for G yields 



(3.34) G(x,t,C) = I + Y t - 

s=i aj 



a 2 G 1 (x,t,(j) <?2 _ 0-3V2G 1 (x,t,Q 0-20-3 

c-o c+o 



Eliminating Aj from (I3.29j) and (I3.32j) . and using Lemma [3.11 together with the 
fact that if)j'_ 1 Dj-.i • • ■ D x = ip^ 1 in the resulting equation, we find that 

(3.35) ^ G ' 1 ( x ^^-(^f^ fa ^)^\x^Q 

for some functions Tj{x, t),Sj(x,t), j = 1, . . . , N. Evaluating equation (13.341) at 
( = (k, and using (I3.35P we find 

(3.36) 

G(li( , &) = I+ gI[ 5 i_g(M)) {bj b -^ (xMj) 

1 frj(x,t) 



-0-3 



Ck + Cj 

But the symmetry (13. 4p together with (13.351) shows that 



%,t)) 6 i 1 )^- 1 (^*, G>3 



(3.37) G(x, t, CO = ^2 (^|) S* 1 )^ 1 ^*,^, 
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so that multiplying equation (13.36!) by C2^o{x, t, ( k ) (b k 1 —bk) T from the right, 
the left-hand side disappears and we are left with 

(3.38) 



=a2ipo{x,t,(k) 
1 

-cr 3 



N 



3=1 



.Cfe - 



r*(a;,t) 



C* + C; W x >*) 



cr 2 ^oOM,Cfc) 



These are 2N algebraic scalar equations for the 2N unknowns {rj(x, t), sj(x, t)}± . 
Given the seed eigenfunction i[) , the solution of (I3.38|) expresses the r/s and the 
Sj's in terms of {Q, bj}f . Thus, by ( I3.32p and (I3.35p . the A/s are determined in 
terms of the 2N complex constants {Q, bj}± . This yields the following result. 

Proposition 3.2. Let N be a positive integer and let {Cj'j &j'}>=i be nonzero 
complex constants such that Q ^ ( k for j ^ k. Assume that u (x,t) satisfies 
equation 112. 1\) and let ipo(x,t,() be an associated eigenfunction obeying the sym- 
metries ( 2.4) ■ Then the following function u(x,t) is also a solution of equation 



N 



(3.39) 



u(x,t) = 2^2(A j (x,t)) 12 + u (x,t), 
i=i 

where the 2 x 2-matrix valued function Aj(x,t) is given by 



(3.40) 



M*,t) = ^ (iffcfl ft b f) ^frAQ 



the constants {cij}^ are given by $3.33\) . and {rj(x, t), Sj(x, £)}fLi are determined 
by the linear algebraic system of equations k3. 38\) . 



4. Solitons 

In this section we apply the dressing method described in the previous section 
to the particular seed solution uq = with corresponding eigenfunction ip = 
exp(—i8a 3 ). The addition of 2N poles to i/j yields the iV-soliton solution of (12. ip . 
Since in this case the solution of the algebraic system (I3.38P can be expressed in 
a simple form, we find simple explicit formulas for the solitons. 

Recall that 6 = 6(x,t,Q was defined in (12. 5ft . Substituting v,o = and 
ip® = exp(—i9a 3 ) into the equations (I3.39P and (13.401) . we find 

N 

(4.1) «(x,t) = 2^-^-, 

where fj := fj(x,t) = bje %e<yX ' t,c - j \ j = 1,...,N. On the other hand, for this 
choice of if) , the top entry of equation (I3.38j) yields 
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where the entries of the N x iV-matrix K = K(x, t) are defined by 
(4.3) K jk = -^(QffE + CO, j,k = l,...,N. 

Combining (14. ip and (14.21) . we find the following expression for the iV-soliton 
solution. 

Proposition 4.1. The N-soliton solution UN{x,t) of the generalized NLS equa- 
tion Ii2.1\) is given explicitly by 



N 

(4-4) u N (x,t) = £ j2(Jr%, 

k,j=i 

where the N x N -matrix K = K(x,t) is defined by the coefficients {fj}i 

are given by 

(4.5) fi := f 3 (x, t) = b, exp ^x + i - t\ , j = l,...,N, 

and the solution depends on the 2N complex parameters {bj,(j}i- 

5. DNLS solitons 
The derivative nonlinear Schrodinger equation ( 11.2jl admits the Lax pair [I] 



(5.1 



i) t + 2i^a 3 ^j = (jrQ - £Q 2 a 3 - \{iQ x a z - Q 3 )) V, 

where the 2 x 2-matrix valued function Q(x,t) is given by 

g N 



(5.2) « - . -j 

If we identify Q with U x , then the x-part of this Lax pair coincides with the 
x-part of the Lax pair (I2.2p . A similar analysis that led to expression (14. 4 p for 
un applied to equations (13 . 1 0[) and (I3.18P yields (cf. [3]) 



(5.3) u Nx 

where K and the f/s are defined as in Proposition 14.11 In view of the identi- 
fication of Q with U x , it is natural to expect the iV-soliton solution of the 
DNLS equation to be given by an expression similar to the right-hand side of 
( 15. 3p . In fact, it was shown in [3] that qN(x, t) is given by the right-hand side of 
(15. 3p with fj(x,t) replaced by 

(5.4) fj(x,t) = b j e i ^ x+2C i At \ 

Now equations (I4.4p and (I5.3P imply that the right-hand side of (15.31) is the 
x-derivative of the right-hand side of (14. 4p when fj is given by (14. 5p . Since all 
the time-dependence in these equations lies in the exponents of fj and fj and is 
irrelevant for the differentiation with respect to x, this property holds also when 
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the /j's are given by (j5.4p . This yields the following simple expressions for the 
DNLS iV-solitons. 

Proposition 5.1. The N-soliton solution q^(x,t) of the DNLS equation ( f 
is given explicitly by 



It can be checked explicitly for small values of N that the solitons and q^ 
given by ( I4.4j) and (I5.5P indeed satisfy the generalized NLS equation (12. ip and 
the DNLS equation (jl.2p . respectively. 



We have implemented the dressing method to equation (12. ip . which is equiv- 
alent to the generalized NLS equation ( II. ip . This provides a way of generating 
new solutions from already known ones. In particular, starting with the trivial 
solution uq = 0, we arrived at the explicit expression (14. 4p for the iV-soliton 
solution. Since (II. ip is related by a gauge transformation to a member of the 
DNLS hierarchy, the construction bears similarities with the implementation of 
the dressing method to DNLS [3J (see also [IH])- The presence of additional sin- 
gularities in the t-part of the Lax pair makes the argument for ( 12. ip somewhat 
more complicated. 

The link between equations (II. ip and (II. 2p allowed us to find a simplifica- 
tion (see equation ( 15. 5p ) of the formulas for the DNLS iV-solitons presented in 
[3]. The fact that the right-hand side of (15. 5p can be expressed as a partial 
x-derivative is related to the fact that the DNLS hierarchy can be reformulated 
in terms of the potential u = d~ l q. Consequently, the solitons of members of 
this hierarchy can be expressed as x- derivatives of elementary functions cf. [9]. 

Acknowledgement The author thanks the referees for valuable suggestions, A. 
S. Fokas for helpful discussions, and T. Tsuchida for bringing the reference [9] 
to his attention. 
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